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Current Acceleration from Dilaton and Stringy Cold Dark Matter 
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We argue that string theory has all the ingredients to provide us with candidates for the cold dark 
matter and explain the current acceleration of our Universe. In any generic string compactification 
the dilaton plays an important role as it couples to the Standard Model and other heavy non- 
relativistic degrees of freedom such as the string winding modes and wrapped branes, we collectively 
call them stringy cold dark matter. These couplings are non-universal which results in an interesting 
dynamics for a rolling dilaton. Initially, its potential can track radiation and matter while beginning 
to dominate the dynamics recently, triggering a phase of acceleration. This scenario can be realized 
as long as the dilaton also couples strongly to some heavy modes. We furnish examples of such 
modes. We provide analytical and numerical results and compare them with the current supernovae 
result. This favors certain stringy candidates. 

PACS numbers: 98.80.Cq. 



I. INTRODUCTION AND SUMMARY 

If string theory is the right paradigm for quan- 
tum gravity then it must provide us with all the in- 
gredients to make the currently observed Universe. 
According to the most recent measurements 
the universe is to a good approximation spatially 
flat, and the current energy budget of the Universe 
consists of about 30% non-relativistic matter, most 
of it being weakly interacting and known as the 
cold dark matter (CDM), and about 70% is in the 
form of an unknown source of dark energy which 
is responsible for the current acceleration. Unless 
a dominant fraction of them is hidden in dark as- 
trophysical objects. Standard Model (SM) baryons 
can only make up a small fraction of the total en- 
ergy density. 

Obtaining the SM sector from string theory has 
been a hot pursuit, see for a review |2]. However in 
making the present Universe it is equally challeng- 
ing to seek a candidate for the cold dark matter 
which cannot be otherwise explained within the 
SM Finally and most importantly, one also 

has to obtain the vacuum/dark energy which seems 
responsible for the current acceleration P, 0| . 

Now it is well known that any string and/or su- 
pergravity compactification to four space-time di- 
mensions usually Q leads to a runaway poten- 
tial for the dilaton (see |^ for efforts to stabilize the 
dilaton) . Let us study the rolling dilaton dynamics 
in the light of dilaton's coupling to various matter 



fields. It is important to realize that the dilaton 
does not have a universal coupling to all the mat- 
ter degrees of freedom , and it is this fact which 
we are going to exploit in our paper j4^ . In or- 
der to have a holistic cosmological evolution in our 
model, we would require the following ingredients: 

1. The CDM candidate should be treated on 
an equal footing as the SM baryons in order 
to explain why flB/^CDM ~ There- 
fore, if we assume that the dilaton couples to 
both the sectors, then it should have a simi- 
lar coupling. In summary, both the SM and 
the CDM component should couple weakly 
(or be uncoupled) to the dilaton, in order 
to comply with various empirical constraints 
from fifth force experiments, from observa- 
tional bounds on the variation of the fine 
structure constant and of Newton's constant, 
etc. 

2. In order to explain the current accelera- 
tion, we further require another CDM sec- 
tor which is strongly coupled to the dilaton, 
and which we coin as SCDM. We will see 
that, by virtue of its coupling to the dila- 
ton, the SCDM redshifts much more slowly 
than CDM and therefore may dominate the 
matter energy density at late times even if 
it was subdominant at early times. More- 
over, if the coupling is strong enough, it slows 
down the rolling of the dilaton, causing the 
energy-momentum tensor of the dilaton to 
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lead to an accelerated expansion of the uni- 
verse. 

We will provide examples of both standard CDM 
and SCDM in the context of string theory. At 
first instance it might look ad hoc to introduce a 
new form of dark matter which is strongly coupled 
to the dilaton, but as we shall see such states are 
already present in string theory. 

Moreover, our approach to explaining dark en- 
ergy has an obvious advantage: we do not have 
to invoke a new scale such as 10"^'^ eV for the 
mass of a new scalar field introduced in order to 
explain the origin of dark energy as in the case 
in models with a quintessence field 0, 0, 0| 
(for alternative explanations, see The natu- 

ral scale in our approach is solely governed by the 
string scale. In this respect our model is minimal 
and facilitates making connections to string the- 
ory. The reason why acceleration starts at a much 
lower scale than the Planck scale, as we argue, can 
be traced back to having small initial cosmological 
energy densities of SCDM string states in the early 
(post-inflationary) universe. 

We begin our discussion with the string theoretic 
motivation for CDM and SCDM (Section 2). We 
explain how the string and brane degrees of free- 
dom we make use of (which are the same ones stud- 
ied in brane gas cosmology see also 
couple to the dilaton, and give examples of can- 
didates for SCDM which couple to the dilaton in 
a way which can give rise to acceleration. In Sec- 
tions 3 and 4, we provide analytical and numerical 
descriptions, respectively, of the evolution of the 
four distinct components in our cosmology, namely 
radiation, CDM, SCDM and the dilaton. We also 
constrain the couplings by fitting to the supernovae 
data. In Section 5, we re- address the cosmic accel- 
eration problem within the string theoretic frame- 
work. We conclude by summarizing and pointing 
to future research challenges. 



II. MOTIVATION FROM STRING 
THEORY 

A. Dilaton Gravity 

The main aim of this section will be to find two 
important couplings, first the dilaton self coupling 



which will determine the slope of the runaway po- 
tential, and second the coupling of the dilaton to 
the stringy degrees of freedom, allowing us to iden- 
tify stringy candidates for CDM and SCDM. 

Let us start with 10-dimensional type II string 
theory whose low energy bosonic action is given by 
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where hatted quantities denote the full higher 
{D = 10) dimensional objects, (p is the dilaton, 
i?3 is the field strength of the antisymmetric ten- 
sor field, and the field strengths F„'s are n forms 
(where n — 2,4 in type IIA theory while n=l,3,5 
in type IIB theory). To be general, we have also 
added a potential for the dilaton. Such a potential 
could result from quantum corrections 0, 0] . 

We will assume that all shape and volume mod- 
uli of the internal manifold are stabilized. This 
could occur by making use of the GKP scenario 
to stabilize the complex structure (shape) moduli 
with the help of various form fields, and then in- 
voking the KKLT ideas Ts'l to stabilize the volume 
modulus with the help of non-pcrturbative pro- 
cesses. Alternatively, but only in heterotic string 
theory, we could use ideas from brane gas cosmol- 
ogy [il in to stabilize both volume |19ll20ll2H 
and shape |22 | moduli by means of string states 
which are massless at the self-dual radius (see also 
for a study of volume stabilization in the string 
frame which is applicable both to type II and het- 
erotic string theories). We assume that these mod- 
uli do not play any substantial role in the late 
time dynamics (this assumption is supported by 
the analyses of 20, 21j), and henceforth we will 
ignore them. 

Performing the well-known conformal transfor- 
mation (relating the string-frame metric on the left 
hand side to a function of the dilaton multiplying 
the Einstein frame metric on the right hand side) 



„40/(D+2) 



(2) 



(where D is the number of compactified spatial 
dimensions), we recover from an action in the 
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10-dimensional Einstein frame 
1 



S = 
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(3) 
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where the new potential is related to the original 
one by rescaling: 



V{<t>) 



(4) 



We take the full higher dimensional metric to be: 

(5) 
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where r is the radius (volume) modulus of the com- 
pact space and where we used the symbol "o" to 
indicate extra dimensional quantities. 

Then, after integrating over the extra dimen- 
sions, one straight-forwardly obtains: 



r 



where we have defined 



V-2y(0)] , 



m: 
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as the four space-time dimensional Planck mass, 
and we have also rescaled 6 



(8) 



in order that it has a canonical kinetic term. 

As in 01, we now specialize to the case when 
V{(j)) is an exponential potential 



y(0) = . 



(9) 



Indeed, such exponential potentials are common 
in lower-dimensional supergravity models and can 
also arise in string theory from quantum correc- 
tions (both perturbatively, for example involving 
string loops 0, or non-perturbatively, like in 
gaugino condensation Although we keep /3 as 
a free parameter in our model, ideally one should 
be able to derive its value from the origin of the 
potential, and we discuss an example in Section 5. 



B. Wrapped Branes as candidates for CDM 
and SCDM 

Branes and strings can wrap certain cycles in 
the internal dimensions. In general, they are heavy 
objects with mass scales of the order of the string 
scale (an exception are the special modes which 
are massless at the self-dual radius - however, such 
modes exist in heterotic but not in Type II string 
theory). The dilaton in general couples to these ob- 
jects, and here we determine these couplin gs. A s in 
the usual brane-gas cosmology paradigm [l^ ITU, 
we assume that branes have annihilated in the 
three non-compact directions, allowing them to ex- 
pand, and leaving us only with branes wrapping 
some or all of the six internal extra dimensions. 
Hence we only need to compute the couplings for 
p-branes with p — 2 ... 6, and for the string wind- 
ing and momentum modes. 

The brane couplings can be derived [251 
starting from the DBI action, which in the string 
frame reads 



Sdbi = 



(10) 



where 7 is the metric on the brane world sheet 
labeled by a. One has to now make conformal 
transformations (|2l followed by dilaton rescaling 
JHl to obtain the brane coupling /ip in the Einstein 
frame: 



S 



DBI 



with 



2/ip 



n/8 



(11) 



(12) 



The wrapped branes of course appear to us, from 
the large four space-time dimensional point of 
view, as point-like objects. Thus, their number 
density redshifts as a~^, like non-relativistic dust. 
In the brane-gas approximation the four dimen- 
sional action for wrapped p-branes are therefore 
given by (see e.g. pol |') 




(13) 
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where pj, corresponds to the observed four- 
dimensional energy densities of the wrapped 
branes which appear in the Einstein equations |43 ■ 
Note that, because of the couphng, if fip is pos- 
itive the energy density Pp redshifts slower than 
if the dilaton increases, and this is why these 
wrapped branes can act as candidates for SCDM. 

By an analogous prescription we can also ob- 
tain the four dimensional action for string winding 
modes [s^, and momentum modes starting from 
the Polyakov action '2^. For the winding modes 
we find 

2.. = ^, (14) 

while the momentum modes do not couple to the 
dilaton at all (27i |: 

Mmom = 0. (15) 

One immediately realizes that, since the momen- 
tum string modes do not couple to the dilaton, 
their energy density redshifts in the same way as 
baryons, i.e. as a~^, and therefore they can be can- 
didates for CDM. On the other hand, the coupled 
winding states can form the SCDM. However, since 
we want SCDM to become important (comparable 
to CDM) only at the present cosmological epoch, 
this means that in the early universe the energy 
density of the winding states must have been much 
smaller than that of the momentum modes. Here, 
we present a qualitative argument based on the 
masses of these objects as to why we expect this 
to be the case. We will address the issue more 
quantitatively in Section 5.2. 

From Eq. (|12|) we can directly read off the masses 
of the p-branes as 

TOp~Mf+VP, (16) 

where r is the radion field in the Einstein frame. 
We assume that the radion is stabilized. Similarly, 
from the Polyakov action the masses for the string 
winding and momentum modes are given by 

mi - M^r , (17) 

and 

1 / ^ 

"Imom ~ - • (18) 



For the purpose of illustration, here we have as- 
sumed that, initially, 4> evolves from close to zero, 
but similar arguements can be made for more gen- 
eral initial conditions. 

If the radion is stabilized at a scale slightly lower 
than the Planck scale, say ~ Mcut, then be- 
cause of the r dependence of the different masses, 
the momentum modes turn out to be the lightest, 
while all the winding string modes and branes are 
much heavier. Thus, if there was a period of in- 
flation of the three large spatial dimensions, and 
if the temperature after reheating is of the GUT 
scale, then the reheating process at the end of infla- 
tion can either perturbatively |2^ or via para- 
metric resonance jsO, 31] produce string momen- 
tum modes 0. The thermal production of the 
winding states (which are very heavy compared to 
the scale of inflation) is exponentially suppressed. 
There are, however, some mechanisms by which 
such heavy modes could be produced [s^, Is^ - 
although their number densities would be much 
smaller than that of the momentum modes. Thus, 
in the context of inflationary cosmology, this may 
yield one way to explain the difference in the initial 
post-inflationary energy densities of the momen- 
tum string modes (CDM) and the winding modes 
(SCDMs). It is possible, however, that the wind- 
ing mode production processes after inflation are 
too weak. In this case, one must invoke the primor- 
dial pre- inflationary density of SCDM states to de- 
termine the late-time number density: in Section 
5 we argue how a period of inflation of appropri- 
ate length can dilute primordial energy density of 
SCDMs to make it just right to address the coin- 
cidence problem. 

If both baryons and CDM momentum states are 
produced during reheating after inflation, then it 
should be expected that their densities are compa- 
rable, as long as the masses of both types of states 
are small compared to the typical energy scale 
during reheating. If the initial post-inflationary 
baryon density is not changed substantially by 
later stages of baryon production, the number den- 
sities of baryons and CDM particles are expected 
to be comparable. 

Let us then assume that after inflation we have 
very small energy densities of the different wrapped 
brane/string gases. Now observe that because of 
the difference in the coupling exponents /ip's, all 
the species redshift differently. As cj) is increasing, 
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it is clear that very soon the largest energy density 
among these different species would correspond to 
the one with the largest positive fip. From (|12|l 
and ((TH) . one finds that this corresponds to the 
6-branes, with 



tion reads: 



Me 



2^8 ■ 



(19) 



Thus, according to our scenario, at late times the 
dark sector of our universe will consist of three 
main components: (a) the dilaton behaving as 
dark energy, (b) string momentum modes play- 
ing the role of ordinary CDM and (c) 6-branes, 
which wrap the six extra dimensions, playing the 
role of SCDM. As we shall see in the next section, 
by virtue of their coupling (|19|l to the dilaton, 6- 
branes can trigger a late phase of acceleration. 



III. COSMOLOGICAL MODEL 



Acceleration due to Coupling: Basic 
Mechanism 



We begin by analyzing a simple scalar-tensor 
theory action for dilaton gravity, where the dila- 
ton has a run-away potential: 




(20) 



In principle, (j) can be identified with any suitable 
moduli field as exponential potentials are quite 
common for other kinds of moduli fields as well, 
for example, in flux compactifications .3,5.] ■ How- 
ever, for the purpose of this paper we focus on the 
dilaton. As in conventional cosmology, gravity is 
sourced by the matter-radiation stress energy ten- 
sor which can be derived from an action of the 
form 



(21) 



where the index i represents the different compo- 
nents of matter/radiation and pi is the observed 
energy density that gravitates. The Hubble equa- 



(22) 

Note there is a departure from the conventional 
cosmology. We allow different matter components 
to couple to (/>: 



(23) 



where we further assume that the ^-independent 
"bare density" pi behaves like a perfect fluid i.e., 
it obeys the continuity equation: 



3i7(p, +p,) =0, 



(24) 



with an equation of state, pi = LOiPi. 

If a scalar field (j) is to act as a quintessence 
type field, then it must be virtually massless, i.e. 

^ 10~^'^ei;. If (j) is coupled to SM parti- 
cles, it mediates a fifth force. Fifth force ex- 
periments/observations constrain such couplings 
of 4) to SM degrees of freedom to be very small, 
MsM < 10^"* — 10^''' (see e.g. l36j). For a dilatonic 
0, a "least coupling" principle was proposed by 
Damour and Polyakov 9] which can explain how 
its couplings to SM particles can become increas- 
ingly small as rolls towards the strong coupling 
regime. We note in passing that observations and 
experiments on variations of physical constants like 
the fine structure and Newton's constants also give 
us similar bounds on p^n- Therefore in the follow- 
ing we will assume that radiation (r), baryons (b) 
and CDM (string momentum modes) have tiny/no 
couplings to (j). 

The dynamics of 4> depends on the exponent /3, 
and the coupling of to various matter-radiation 
components. The equation of motion for is given 
by 



with 



Kff(0) = Foe-2/9^ + 



E 



(25) 



(26) 



where the sum runs over all matter states. For 
simplicity, we have assumed that the gauge bosons 
do not contribute to an effective potential for (p. 
This is strictly true if on average the electric field 
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and magnetic field vanishes. There could be other 
forms of radiative matter which could induce an 
effective coupling, as discussed for example in |3^ . 
Weakly coupled matter (stringy candidates for 
CDM and baryons) are denoted by the subscript 
"m", and the SCDM by the subscript "sc". 

It was argued in [2J| that, when the Universe 
evolves adiabatically, the dilaton (f) tracks the min- 
imum of the effective potential which evolves in 
time with the redshifting of the p^'s. How fast or 
slow the dilaton is evolving mostly depends on the 
equation of state parameter Ui, and on the cou- 
pling exponent /ij of the dominating matter com- 
ponent. Thus, at earlier times when radiation or 
weakly coupled matter is dominant, (f) rolls fast and 
tracks the relevant matter/radiation component, 
remaining subdominant. However, once SCDM be- 
comes important, the dilaton slows down, starts to 
dominate the energy density, and our universe en- 
ters a phase of acceleration. Let us now see this in 
more detail. 



B. Exact Two-fluid Analysis 

We will first work within the approximation that 
there are only two dominating components in the 
Universe, the dark energy determined by (j) and a 
single matter/radiative component. To keep things 
general, we do not specify which form of matter it 
is. In this approximation, the Einstein constraint 
equation l|22|) determining the Hubble expansion 
rate reads 

^ (27) 



3 V ' 2' 

From Eq. (|25() . it follows that the evolution equa- 
tion for (j) is given by 

4> + iH4> = 2{[3VQe-^^'^ ~ Aie^^'^p) ■ (28) 
Solving the continuity equation H24|) . we obtain, 

-3(1+1.^) 

(29) 



P = Po — 

where po is the present day value of the energy 
density. 

Performing a change of variables r\ = ln(a(t)) 
leaves us with the following set of equations: 



1^ 



2/3|^e-^^ 



where ' = d/drj. The system of equations can be 
solved analytically to obtain the "tracking solu- 
tion" , the solution where 4'{t) tracks the minimum 
of its "effective" potential. The solution is given 
by 



/3(/3 + /i)- 3/4(1 + a>) 



Vo 

p \pip + (3)+ 3/4(1 + u)^ - 3/8(1 + Lu) 



(30) 

p 



where 

and 

a(t) = ao 



1 



P = 



2{f3 + p) 



(31) 



with ujfj, — 



p/13- 



plf3 + l ■ 
(32) 

A few remarks are now in order. First, we ob- 
serve that the energy densities p and p^ track each 
other and redshift with precisely the same equa- 
tion of state parameter oj^ that was derived in the 
adiabatic approximation in p3 |: 



P<t> = 



pin + /3)+ 3/4(1 + uj)^ ~ 3/8(1 + u) . 
Pif3 + p)-S/A{l + uj) 

-3(l+w^) 



-P 



rp^a ' . (33) 
One recovers the value of the tracking ratio 



(34) 



obtained in the adiabatic approximation j24| when 
P-T f3 > C'(l) Bind the terms involving (1 + uj) can 
be ignored. 

Second, we realize that real solutions for only 
exist if the numerator on right hand side of 
is positive (the denominator is positive as long as 
u! > —1/2). In other words, in order for a tracking 
solution to exist we require: 



PiP + p) > 



3{l + to) 



(35) 



Since we want such tracking behaviour to hold 
during the radiation-dominated phase, substitut- 
ing pr — and = 1/3 in Eq. (|35l) . we find: 



/3 > 1. 



(36) 
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We note in passing that, if the condition Eq. (|35|l 
is not satisfied, then we end up in a different non- 
tracking late time attractor solution. In this case 
the field does not follow the minimum but lags 
more and more behind as the universe expands, 
and the energy density becomes dominated by the 
dilaton from the very beginning. From the point of 
view of cosmology this is an uninteresting scenario. 

Finally, we find from Eqs. H32|) that acceleration 
happens when 



-(- 



(37) 



For winding brane SCDM {lOsc — 0), combining 
l)37|l with H36(l . we find that we can account for a 
late acceleration phase provided our SCDM has a 
coupling exponent 



Psc > 



(38) 



Note, in particular, that the exponent which 
we computed for wrapped 6-branes satisfies this 
bound. 



C. Three fluids system 

Our exact analysis of two fluids model suggests 
that, if the dilaton couples to SCDM, then there 
can be an interesting dynamics which triggers a 
phase of acceleration. In order to study more 
precisely how such a transition happens, one has 
to consider a more realistic cosmological model 
where, in addition to the SCDM fluid, we also have 
a normal CDM component that is uncoupled or 
weakly coupled to the dilaton. Such a model has 
the advantage that, at early times, the evolution 
is standard, i.e. H ~ a~^^^, and, at late times, 
as SCDM plus dark energy begin to dominate, we 
can have accelerated expansion. 

The new system of equations (|22I25|) becomes: 



it hardly makes any difference to our analysis. Al- 
though no analytical solution to Eqs. (|39I40|) is 
found, one can obtain some insight into the cos- 
mological evolution resulting from these equations 
by considering various limiting situations. 

First transition: First, consider a situation 
when the energy density in SCDM is much smaller 
than that in both ordinary matter and in V{(f)), so 
that one can neglect the second terms in Eqs. (|39fl 
and (|40|l . In this case, the equations from the two 
fluid discussion apply, making use of /Lt = and 
w = 0. This is the phase when the dilaton es- 
sentially rolls down freely along its potential V{(t>), 
tracking the energy density of ordinary matter |24| 
as given by l|55|) : 



8/32 



6 



Pn 



(41) 



The scale factor a{t) cx t^^^ as in the usual matter 
dominated era. Note that during this period the 
value of (f) which is the minimum of the "effective 
potential" for (j> (the sum of V{4i) and psc) is red- 
shifting more slowly than the value of (p{t) in this 
phase (this can be seen by going back to H30|l , since 
it is setting ^ = in this equation which gives the 
time dependence of the minimum), and thus (j){t) 
eventually catches up with the minimum. After 
this point, the value of (f)(t) tracks the minimum. 
In this case, we know that the energy density of (fi 
will track SCDM, as given in Eq. H33() . It is clear 
then, that when the right hand sides of Eqs. (|33|l 
and (|69|1 become equal, i.e. when 



4/32 



P' 



and thus 



3/3 



M(4/32-3)' 



(42) 



(43) 



a first transition, namely the transition between 
the period when (/)(t) is moving freely and when it 
I starts tracking the minimum of the effective poten- 

+ Psc + Pm) , (39)tial (set by the interaction of the dilaton potential 

and the contribution from the SCDM term) occurs 
— PscPsc^'^^'"'^ — MmPmS^'^™'^) j^. From this moment on, both psc and pq start 
- PscPsc - PmPm] , (40)*° redshift as Eq. (jSSl, with w -> and ^ 

Second Transition: The second transition |53| 
For simplicity, we assume for the time being that corresponds to when the universe enters a phase of 
pm = 0; we will see later that as long as << (3, acceleration. This happens approximately when 



3Mp2 
= 2[/3y((/.)- 
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Pd = Psc + P4, = (1 + rsc)psc becomes equal to 
the energy density of the ordinary matter compo- 
nents, as can be readily seen from the Hubble rate 
Eq. Thus, we enter a phase of acceleration 

when: 



Ps 



1 



1 



1 



-Pn 



(44) 



where in the last step we again made use of the 
adiabatic approximation. 

First, notice that depending on how large the 
ratio Psc/ P is, one needs a very small fraction of 
SCDM compared to ordinary weakly coupled mat- 
ter to initiate the phase of acceleration. Since 
SCDM behaves like non-relativistic dust {lOsc — 0) 
it would cluster around galaxies and contribute 
to the current dark matter abundance, thereby 
changing the baryon to dark matter ratio in re- 
cent times. Since there is reasonable agreement 
between this ratio as obtained from the cosmic 
microwave background (CMB) (which measures 
the primordial fluctuations at redshifts of around 
z = 1100) and measurements of rich clusters of 
galaxies (at recent redshifts), it is important that 
this ratio does not change too much in our model. 
This is guaranteed if p^r. is only a small fraction of 
As pointed out in j2J| this scenario of chang- 
ing baryon to dark matter abundance offers the 
intriguing possibility of reconciling a 10 to 20 per- 
cent discrepancy in the estimates of baryonic abun- 
dances coming from CMB and big bang nucleosyn- 
thesis (BEN) measurements |37l| . 

What if SCDM is a radiative fluid?: So 
far, we have assumed that SCDM is a form of dark 
matter (winding states) but is strongly coupled to 
the dilaton. What if it were to behave as a radia- 
tive fluid? For example, in j38j | it was proposed 
how thermal fluctuations of a massless field with 
quartic interactions can act as a radiative fluid 
with precisely the kind of exponential coupling to 
the dilaton that we considered here. In this case 
= 1/3, and we will observe it as a component 
of dark energy. This is because after the first tran- 
sition both SCDM and (/) obey the same equation 
of state as computed in Eq. Ip!!!jl , and cosmolog- 
ically we cannot distinguish between the two. The 
total dark energy would be then given by pd- Obvi- 
ously, in this setup the ratio between dark matter 
and baryons remains constant and we do not en- 
counter any constraints on the ratio Psc/ P on this 



basis. From we find that in order to have 
acceleration the ratio just needs to satisfy 



(45) 



Acceleration epoch: In either case, the main 
constraint comes from trying to match the ob- 
served equation of state, uj^ < —0.76 From 
Eq. (|32|l we find such a bound implies 



for SCDM and 



f3 



> 3 



> 5 



(46) 



(47) 



for radiation. How to get such large value of the ra- 
tio from string theory is discussed in the following 
section. 

For these ratios of Psc/ P one can also estimate 
approximately the onset of the period of accelera- 
tion. It is the time when pd ~ Pm- Now 



1 = 



Pd.acc 




^acc \ 


Pm^acc 




V ao / 


^dQ f 


1 ^ 






l~ ^acc J 





(48) 



Using the known dark matter to dark energy ratio, 
we find 



0.5-1.0, 



(49) 



which is in good agreement with the supernova 
data. 

Finally, we note that most of our analysis of the 
transitions go through identically even if pm is non- 
zero but small. This is because, if pm is tiny it 
changes very little the tracking ratio in the matter 
phase (|69|) which enters in the calculation for the 
epoch of first transition. Also, 



Pn 



-3+0(p™//3) 



Thus, weakly coupled matter redshifts essentially 
the same way as in the standard cosmology and 
therefore the second transition is also unaffected. 

To summarize, we have three distinct cosmolog- 
ical eras depending on which terms are most domi- 
nant on the right hand side of the Hubble equation 
Eq. (1231: 
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1. Radiation era, when psc and pm are negligi- 
ble while p0 tracks p^i 

2. Matter era, when and are negligible 
while p0 tracks and finally, 

3. an Acceleration phase, where stringy compo- 
nents, viz. starts to dominate over pm. 

As far as the parameters of our model are con- 
cerned, first we find that in order to have a late 
time acceleration phase we need the ratio /Zsc//3 
to satisfy conditions like Eq. H70() . On the other 
hand, to have an earlier non-accelerating phase 
fir , pm has to be small which is consistent with 
constraints coming from Fifth force experiments, 
and observations of variation of Gjq and a / which 
lead to the condition psM < 10"^. Finally, the 
requirement that (f> should track radiation in the 
early universe tells us that (3 > 1. In spite of these 
constraints, there is a huge parameter space where 
we seem to agree with the cosmological and parti- 
cle physics observations, which is encouraging. Let 
us therefore study this model in more detail using 
numerical techniques and see whether the allowed 
parameter space can be made more precise. 



IV. NUMERICAL WORK 

A. Evolution and Abundances 

In order to verify our analytical approximations 
and see the precise evolution of the system, we 
solve the three fluid system numerically. We start 
deep in the matter dominated era and choose ini- 
tial conditions such that the field initially follows 
the approximate early time solution. As an exam- 
ple, the evolution of the field for the set of param- 
eter values (3 = 6, p= 10, flrh = 0.27, = 0.01, 
and Vo fixed by requiring flatness is shown in Fig.Q] 
The plot depicts the field value as a function of 
time (measured in terms of the logarithm of the 
scale factor). The solid red curve represents the 
full numerically determined evolution. In the same 
figure we also show the analytical approximations 
given by Eqs. (|39I4U|) in a universe that is radiation 
dominated, i.e. exp(0) = {SP^Vo/SprY^^^ (dot- 
dashed purple curve), dominated by ordinary mat- 
ter, i.e. (exp(0) = (8/32Vb/3/5„)i/2/3 ^^^^ dashed 



0.3 p 
0.2 - 
0.1 - 




log(a) 

FIG. 1: The evolution of as a function of the log- 
arithm of the cosmological scale factor for the full nu- 
merical solution (solid red) and for the approximations 
to the solutions corresponding to the evolution at late 
(blue dashed) times, early times (green long dashed) 
time and in the initial radiation-dominated (purple 
dot-dashed) phase. 

1 I 1 1 1 1 I 1 



0.8 - 
0.6 - 




-6 -5 -4 -3 -2 -1 
log(a) 

FIG. 2: The relative energy densities in the (j> field 
(red solid), in coupled CDM (SCDM)(blue dashed), 
uncoupled CDM (green long dashed) and in radiation 
(purple dot dashed). 



green curve), or described by the late time solu- 
tion (dashed blue curve) when the field follows 
the effective minimum of the potential Eq. 133|l . 
From the figure one can see how initially the field 
follows the radiation dominated solution until it 
crosses the matter dominated solution, after which 
the evolution of the fields is well described by it. 
At late times, the field follows the solution of the 
evolution of the minimum of the effective potential. 
Note also how the crossover from one solution to 
another is very rapid and smooth. 
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25 
20 



15 - 



5 - 



1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 
P 

FIG. 3: Contour plot of Ax^ {^a = 0.27, fi™ = 
0.01). Red (solid) contour corresponds —4, green (long- 
dashed) to —2 and blue (sliort-daslied)to 0. 




(6CDM 
Eds 
ACDM 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

z 

FIG. 4: Luminosity vs. distance with respect to the 
empty universe for (/3 = 1.3, fi — 8) (green, long- 
dashed), ACDM model (purple, dot-dashed) and EdS 
(blue, short-dashed) along with the SNIa data. 



The relative contributions fli of various matter 
constituents (labeled by i) to the critical density 
of a spatially flat universe are shown in Fig. |21 It 
is clear how at late times the coupled dark mat- 
ter component starts dominating along with the (j) 
field. 



B. Observational Constraints: Supernova 
type la data 

In order to see whether the model presented here 
is compatible with present observational data, we 
compare it with the recent expansion history as 



probed by the supernova type la (SNIa) data. Here 
we use the sample of 194 SNIa presented in 39] . In 
order to compare our model with the ACDMmodel, 
we compute the for each set of parameters and 
compare it with the value for the corresponding 
flat ACDMvalue, which for the same set of super- 
novae is Xacdm — 201.6 {rim = 0.27, JIa = 0.73). 
The contribution from radiation is neglected for 
the SNIa fit, and the universe is assumed to be 
flat. The Hubble parameter is marginalized over 
in all of the fits (by maximizing the likelihood) . 

In Fig. Owe show a contour plot of coefficients 
{fi, (3) corresponding to Ax^ = X^-x\cdm^ where 
we have fixed = 0.27, Vise = 0.01. The re- 
gion above the blue (short-dashed) curve shows 
where the is smaller than Xacdm ^ indicating 
a nominally better fit to the data. The green and 
red curves indicate contours where Ax^ is —2 and 
—4, respectively. The reduced x^ the ACDM 
model is 201.6/(194 - 1) « 1.045, so an equally 
good reduced is achieved with the model con- 
sidered here (with two extra parameters), when 
Ax^ « -2.1. 

The dependence on the coupled matter parame- 
ter flsc is very weak, whereas the uncoupled com- 
ponent plays a significant role. For example, with 
large values, flm > 0.4, the contours shift such that 
all of the considered {/3, fi) parameter space gives 
a larger x^ than the ACDMvalue. Small values, 
rtm < 0.2, lead to the same conclusion. Here our 
main purpose is to demonstrate the properties of 
the model considered here and to show how it can 
fit the current SNIa data very well. A more com- 
plete parameter scanning is left for future study. 

In order to understand qualitatively why Ax^ is 
negative for a range of parameter values, we show 
the luminosity distance plot as a function of red- 
shift with respect to the empty universe, AM, for 
a set of parameters {/3 — 1.3, = 8), for which 
- Xacdm = -4-2, in Fig. H From the fig- 
ure we see how for these parameters the luminos- 
ity distance at low redshifts closely follows that of 
the ACDM but at later times curves down more 
quickly, fitting the high redshift data points bet- 
ter. 
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V. STRINGY PHYSICS AND COSMIC 
ACCELERATION: 

Although we have seen that stringy SCDM can 
pave the way for a late time acceleration phase, 
challenges remain to connect string theory with 
our coupled quintessence model in a quantitative 
way. Firstly, our numerical analysis involving su- 
pernovae data suggests that we require /x//3 > 2 
and P > 1.4 for a good fit. The parameters are fur- 
ther pushed up if one also considers bounds coming 
from BBN: During the radiation era = 0) we 
have 



6(/32 



:Pr ■ 



(50) 



BBN considerations provide an upper bound on 
the dark energy abundance, fi^ < 0.21 during the 
radiation era. From Eq. H5()|l we see that for this 
to be satisfied one needs /? > 2. 

Secondly, it is clear that the success of our cos- 
mological scenario relies on having a very small 
abundance of CDM and SCDM in the early uni- 
verse, and it would not be much progress if we sim- 
ply shift the usual problem of quintessence mod- 
els, namely the problem of having to invoke a new 
and extremely small mass hierarchy (namely the ev 
mass scale), to a problem of an unnaturally small 
initial abundance of the SCDM particles. Thus in 
order to truly address the cosmic acceleration issue 
one needs to explain quantitatively why we expect 
such small abundances. In this section we address 
both of these issues, starting with the first. 



A. Quantum Stringy Corrections 

One may argue that we can obtain the required 
values of /x, f3 if stringy SCDM is coupled to sev- 
eral scalar fields This is because, depending 
upon which linear combinations are stabilized, the 
effective for (j) (the linear combination which is 
rolling), can range from to y^J2i M/i where ^/'s 
are the coupling exponents to the different scalar 
fields, 0/'s. Analogous arguments hold for /3. How- 
ever, here rather than invoking many scalar fields, 
we furnish an example showing how quantum cor- 
rections involving stringy loops may be able to pro- 
vide us with a significantly large value of fi. 



In the string frame quantum string loop correc- 
tions will in general modify the couplings of the 
dilaton to the different fields: 



S 



d'^'x [C{^)R + K{(j)){d<j)Y - V{(j})] , (51) 



where the functions C {(p) , K {(/)) and V{(t)) are ex- 
pected to have a Taylor series expansion in terms 
of the inverse string coupling constant g^^ = e~'^ 



c(0) = ^"^""^ 



(52) 



n=0 



and similarly for K{4>) and V{4)). We recover the 
classical action ^ when all co-efficients, except 
C2 = 1 and K2 = 4, vanish. Now, let us sup- 
pose that the first non-zero coefficients of C((/)) 
and K{(t)) occur at the n'th level, while Vb = 
[ssj ]. Then, as — > 00, the action is given by 



S: 



d'\ [e-"^i? + e-"^Knidcb)'' - Vie 



n4> 1 



(53) 

To go to the Einstein frame we perform a conformal 
transformation 



9mn 

followed by rescalings 



\ 



D +2 



Kn(p = kcj) , 



(54) 



(55) 



to make the kinetic terms for the dilaton canonical. 
After dimensional reduction to 4 dimensions, the 
action reads 



5* = 



1 



16nG 
where 



dx'^V^[R-d,n(l>d"'(p~Vie'^'^'f'] , 

(56) 



2/3 



D +2-2n 
k{D +2) 



(57) 



Next, we look at the DBI action for £)-branes 
and perform the conformal rescalings Q to obtain 
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the effective 4 dimensional action for the gas of 
6-branes: 



(58) 



ao 




with 



D (n - 1) + 71 - 2 



(59) 



k{D +2) 

Note that we want /i//? to be positive. Substituting 

o 

D= 6 one finds that this is possible only for n = 
2,3,4: 



For n = 2, 



/3 



1 



Forn = 3, ^ = — ^ 3.13 

p 4 



For n = 4, — = 10 

P 

From Eq. (|55|l it is also clear that, depending on 
the value of Kn, k could be small, hence making 
Hscfi large. 



B. Primordial Inflation and the Scale of 
Current Acceleration 

It is clear that a resolution of the dark energy 
problem should not only involve (a) a tracking 
mechanism which explains why the dark energy 
density has always been close to the matter den- 
sity, and not just today (and this is most certainly 
the case in our model), but also (b) a parameter 
(since we do not want to introduce a hierarchical 
ev mass scale) which governs when the dark en- 
ergy comes out of the tracking phase and starts to 
dominate the universe. Moreover, this parameter 
should be such that its expected value can nat- 
urally explain why we are entering this phase of 
acceleration so late in the day. 

Now, in our model, acceleration commences ap- 
proximately when SCDM becomes comparable to 



ordinary CDM. Since SCDM redshifts more slowly 
than CDM, this means that in the early universe 
(say just after inflation) the energy density in 
SCDM was much smaller than that of CDM. This 
could happen, for example, if the SCDM is pro- 
duced very scantily during the reheating process, 
as argued in Section II. More importantly, one re- 
alizes that SCDM candidates are stringy winding 
states in nature. One expects these winding states 
to be present in the very early pre-inflationary uni- 
verse with approximately string scale energy den- 
sities. However, as our universe expands exponen- 
tially during inflation, this gas of branes is going 
to become very dilute. Since the slope of the po- 
tential decreases exponentially as </> increases, we 
expect the dilaton (p to be effectively frozen dur- 
ing the phase of inflation. As a result, the SCDM 
energy density will redshift as . 

Thus, at the end of inflation, while the energy 
density in ordinary matter would be given by the 
reheat temperature, the energy density in the pri- 
mordial gas of SCDMs would be much smaller due 
to inflationary dilution. In fact, the larger the 
number of e-foldings A/", the greater the dilution, 
and therefore the longer it will take for SCDM to 
catch up with CDM, and the smaller would be the 
scale (= Ma) at which the second phase of accel- 
eration commences. If Ma = MplO~"*, then one 
finds (see Appendix for details) 



R{1 + 4/i//3 - iujsc) + (1.37V - 27)(1 + ljsc) 



4(^1//?- 



(60) 

where AfplO ^ corresponds to the reheat temper- 
ature. 

We are now in a position to estimate the scale 
at which the second acceleration phase starts. As 
a first example, let us choose lOsc = 0, and = 
1/2, the minimal value required for acceleration 
and also approximately what one obtains for the 
6-branes. For the minimal number of e-foldings 
A/" = 60 and assuming reheating to GUT scale tem- 
peratures Mreh, i-e. Mreh ~ 10~^Mp Or i? = 3, one 

30 / 



ev. 



finds from GOJ) that Ma ^ lO'^^Mp = 10" 
We remind the reader that the current energy den- 
sity corresponds to the energy scale 10 ~^ ev. For 
LUsc = 1/3, i.e. radiative SCDM, one needs ^/f3=l 
for acceleration. Then using previous values for J\f 
and R we again find Ma ^ lO^'^^Mp. 

The agreement, however, is not as dramatic as it 
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seems. The more detailed numerical analysis per- 
formed with ujsc = revealed that we get good 
agreement with the supernova data for values of 
^ 2. Demanding A ~ 31, and ~ 2 one 
"predicts" the number of primordial inflation to be 
Af ~ 200 if one demands that the second stage of 
acceleration is beginning at the present time. Es- 
timates with uJsc — 1/3, and ^/f3 ^ 5 (BTjl yield 
a similar estimate for Af, namely Af ^ 250. Of 
course, this number depends on the reheat tem- 
perature, and on the precise pre-inflationary initial 
energy density of SCDM, but these dependences 
are relatively weak as can be seen from (I7U|I . Thus, 
l)70|l should be treated as a relation between A and 
Af. Moreover, since A — A{Af) is a linear rela- 
tionship, there is no fine tuning involved; a small 
variation in Af results in only a small variation in 
A. This is of course to be contrasted with the fine 
tuning problem in ACDM models where the differ- 
ence of two large numbers is supposed to yield a 
very small cosmological constant. Hence, a small 
variation in one of the large numbers causes a huge 
variation in A. 



VI. CONCLUSIONS 

In this paper we have introduced the concept 
of "stringy cold dark matter" (SCDM). From the 
point of view of our four-dimensional space-time, 
SCDM are particles which couple non-trivially 
(with a coupling exp{2fi(f))) to the dilaton. Can- 
didates for SCDM from string theory are branes 
wrapping some or all of the compact spatial di- 
mensions. We have studied the cosmology of these 
modes in the context of a Type HA supergrav- 
ity action to which we have added a run-away 
potential Vbea;p(— 2/30) for the dilaton (j) (as ex- 
plained in the text, such potentials are generally 
believed to be generated when going beyond the 
classical description). We have shown the existence 
of cosmological tracking solutions in the context of 
which the dilaton becomes a candidate for dark 
energy [56[. 

Our model contains radiation, regular dark mat- 
ter (which couples very weakly if at all to the dila- 
ton), SCDM and the dilaton. Assuming that our 
four-dimensional space-time underwent a period of 
primordial inflation, the initial energy density in 
SCDM is expected to be exponentially suppressed 



compared to the density of regular CDM and ra- 
diation (which can be produced during reheating 
after inflation). In this context, the universe is ini- 
tially dominated by radiation, and we show that 
the energy density in the dilaton tracks that of 
radiation until the dilaton gets stuck at the mini- 
mum of its effective potential, a potential formed 
from the dilaton potential (potential with a nega- 
tive slope) and the terms with positive slope com- 
ing from the interactions of the dilaton with the 
SCDM. From then on, the density in SCDM and 
in the dilaton redshift less fast than that of or- 
dinary matter. During this phase, the density in 
the dilaton is greater than that in SCDM. Once 
the dilaton energy begins to dominate, a period 
of acceleration begins, provided fi/ (3 is sufficiently 
large. 

Compared to models of quintessence, an advan- 
tage of our approach is that it does not require the 
introduction of a new mass hierarchy (a new mass 
scale of about leV). In order to explain why the 
late-time acceleration begins at the present time, 
a sufficient suppression of the number density of 
SCDM states is required. Such a suppression is 
naturally generated by a period of cosmological in- 
flation. In fact, the number Af of e-foldings of infla- 
tion required to suppress the SCDM density is (for 
GUT-scale reheating) larger but of the same order 
of magnitude as the minimal value of Af required 
to successful inflation. Thus, in our framework the 
"coincidence mystery" of why dark energy is be- 
ginning to dominate today is tied to the duration 
of the period of inflation. 

For our scenario to work, a sufficiently large 
value of (u//3 is required. Even larger values of this 
ratio are required to be consistent with big bang 
nucleosynthesis. At the level of the classical action, 
it is not possible to obtain such large values with 
SCDM particle couphng to a single scalar mod- 
uli. The cumulative effect of SCDM coupling to 
several different scalar fields may help resolve this 
difficulty, as may non-perturbative effects. Further 
research on this issue is required. 

It would be interesting to study further conse- 
quences of the existence of SCDM in the present 
universe, following the approaches of to study 
cosmological consequences and of to study as- 
trophysical and particle physics aspects. 
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try to estimate when SCDM becomes comparable 
to ordinary matter/radiation. 

After inflation, all the scalar fields are free to 
roll. The energy density of the (j> field will first 
track the radiation and then matter energy density, 
maintaining a constant ratio with them. Thus, in 
the radiation era 



VIII. APPENDIX 

We start by assuming that in the very early uni- 
verse the energy density of the stringy SCDM and 
of is given by the string scale, psc ^ V{Q) ~ 
M^j^. Now, a small hierarchy between the string 
scale and the scale of inflation ensures that Q or 
SCDM play no role in the inflationary mechanism 
and that (j) is effectively frozen during this phase. 
As a result, the SCDM energy density redshifts 
as a-3(i+'^sc) -^ith the rapid expansion of our uni- 
verse. Accordingly, after the end of inflation the 
SCDM energy density will be given by 



PscV-R) — Pie 



(61) 



where Af is the number of e-foldings of inflation, 
and tn stands for the time of the end of inflation 
(the time of reheating), and pi is the energy den- 
sity at the beginning of inflation (when the inflaton 
potential becomes dominant). It is a good assump- 
tion that the SCDM density will not be suppressed 
relative to that of other matter between the initial 
time and the onset of inflation. 

We choose, by convention, (j){tji) = and thus 



-2/30 



, or Vq 



Let us also 



assume, just for the purpose of illustration, that 
reheating is very fast, i.e. that the period of infla- 
tion is immediately followed by a radiation period 
(this is likely to be the case of reheating is driven 
by parametric resonance). In this case, after re- 
heating, the radiation energy density is given by 



Pr{tB) = Pi 



4R 



(62) 



where R sets the scale of inflation (which is 
10~^Mp). It is now clear from (EH) and (EH) that 
after the end of inflation for Af > 60, and Wsc > 0, 
the radiation energy density is a lot larger than the 
SCDM density. This is essentially why it takes a 
long time for to catch up with radiation den- 
sity. As wc have seen, once it does, a second phase 
of accelerated expansion begins. Let us therefore 



(63) 



Hence the radiation and SCDM densities redshift 
differently 



4:tl/f3-3{l+LOsa) 



M|10"^ 



.3(l+w„)A/'"4H / Q \ 



and 



Pr = Mho-'" 



ttE 



(64) 



(65) 



After the radiation-matter equality, when the 
energy density of the universe has fallen to about 
10^^°*Mp, (j) starts to track the matter density in- 
stead 



a 



(66) 



Now, one can obtain the ratio a^^/ a^, where "teg" 
is the time of radiation-matter equality and a^q is 
the value of the scale factor at that time, and then 
substitute it in (|64|l to obtain the energy density 
of SCDM at the equality epoch: 

p.,((te,) = Mpno-i-3(i+--)^-4« 

10(27-fl)(4M//3-3(l+(^,,)) ^ ^gy^ 

The evolution of strongly coupled and ordinary 
matter is then given by 



Psc Pscipeq) ( 



3Ai//3~3(l+aj,,) 



and 



p„ = m;io 



4in-108 



(68) 



(69) 



respectively. 
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Now, as explained before, psc ^ Pm = M\ 
corresponds to the acceleration epoch. If Ma = 
MplO"^, then from ^ and ^ we have 

_ R{\ + 4/i//3 - iLUsc) + (1.3AA ^ 27)(1 + cj,,) 
4(/i//3 - Use) 

(70) 
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